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S Abstract 

, , For Bezier curves, subdivision algorithms create control polygons as piecewise 

linear (PL) approximations that converge in terms of Hausdorff distance. We 
prove that the exterior angles of control polygons under subdivision converge 

^ to at the rate of O(yJ^), where i is the number of subdivisions. This an- 

gular convergence is useful for determining topological features. We use it to 

, show homeomorphism between a Bezier curve and its control polygon under 

subdivision. The constructive geometric proofs yield closed-form formulas to 
compute sufficient numbers of subdivision iterations to obtain small exterior 
angles and achieve homeomorphism. 
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1. Introduction 



5_i A Bezier curve (Definition 3.1) is characterized by an indexed set of 



points, which form a PL approximation of the curve, called a control polygon 



(Definition 3.1). The de Casteljau algorithm [8J is a subdivision algorithm 
associated to Bezier curves which recursively generate control polygons more 
closely approximating the curve under Hausdorff distance [21]. 
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There is contemporary interest [TJ [21 [151 EES] to preserve topological char- 
acteristics such as homeomorphism and ambient isotopy between an initial 
geometric model and its approximation. In particular, the control polygon 
homeomorphic to a simple Bezier curve is also simple, so homeomorphism 
precludes undesired self-intersections, while the control polygon ambient iso- 
topic to a Bezier curve has the same knot type as the Bezier curve. These 
are fundamental for applications in computer graphics, computer animation 
and scientific visualization. 

However, there may be substantial topological differences between Bezier 
curves and their control polygons. First of all, Bezier curves and their control 
polygons are not necessarily homeomorphic. There are examples in the lit- 
erature showing simple Bezier curves with self-intersecting control polygons 
or self-intersecting Bezier curves with simple control polygons [TT1 |23| I2H] . 
Secondly, Bezier curves and their control polygons are not necessarily ambi- 
ent isotopic. There is an example showing an unknotted Bezier curve with a 
knotted control polygon 0, |T5]. An example of a knotted Bezier curve with 
an unknotted control polygon was constructed recently [TT] . 

Control polygons converge under Hausdorff distance does not necessarily 
yield the equivalence relations mentioned above. We find that the angular 
convergence, which to the best of our knowledge has not been previously 
established, is useful for determining homeomorphism and knot type. Here, 
we first show the angular convergence and then prove the homeomorphism 
based on the angular convergence. We consider knot type separately and 
present that in the companion paper [12]. 

Computationally, it is known that the convergence in Hausdorff distance is 
exponential [22] . We show that the angular convergence rate is also exponen- 
tial. Readers will find that the algorithmic efficiency of achieving the equiva- 
lence relations depends on these exponential convergence rates. Furthermore, 
we derive closed-form formulas to compute sufficient numbers of subdivision 
iterations to achieve the desired topological characteristics. These formulas 
rely upon the constructive geometric proofs presented here. 

2. Related Work 

The Angular Convergence we show here is based on a previous result 
showing that the discrete derivatives of the control polygons converge expo- 
nentially to the derivatives of the Bezier curve |20j . 
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The topological characteristics for a Bezier curve and the control polygon 
have been studied in the literature. A previous proof [23J of homeomorphism 
under subdivision relied upon the hodograpfj^] and did not provide the num- 
ber of subdivision iterations. We provide a constructive geometric proof that 
sufficiently many subdivisions will produce a control polygon homeomorphic 
to a given Bezier curve, and a sufficient number of subdivision iterations is 
also established. Topologically reliable approximation in terms of homeomor- 
phism of composite Bezier curves was also studied by the paper [6], which 
used algorithmic techniques that do not completely rely upon de Casteljau 
algorithm, techniques with choosing so called "significant points" first. 

Ambient isotopy is a stronger notion of equivalence than homeomorphism. 
There exist a result [9] establishing this equivalence for geometric models, 
by a certain algorithm distinct from subdivision. Besides, recent papers 
[U [13] present algorithms to compute the isotopic polygonal approximation 
for 2D algebraic curves. Computational techniques for establishing isotopy 
and homotopy have been established regarding algorithms for point-cloud 
data by "distance-like functions" [5]. In the companion paper [TJ], the input 
data there is different of a Bezier curve, with the important additional result 
of the number of subdivision iterations needed to attain an ambient isotopic 
approximation of the original Bezier curve. 

Specifically for 3-D Bezier curves under subdivision, the paper [Tj5] showed 
ambient isotopy only for curves of low degree (less than 4), where a crucial 
unknotting condition was trivial. Deriving a comparable unknotting condi- 
tion for higher degrees entailed significant new arguments that depend on 
our angular convergence. 

A geometric technique we use is constructing a tubular neighborhood 
for a Bezier curve, whose boundary is called a pipe surface. Pipe surfaces 
have been studied since the 19th century [T7j, but the presentation here 
follows a contemporary source [14]. These contemporary authors [H] provide 
Bezier curves considered here. These authors perform a thorough analysis 
and description of the end conditions of open curves. The junction points of 
a Bezier curve are merely a special case of that analysis. 



2 The derivative of a Bezier curve is also expressed as a Bezier curve, known as the 
hodograph [5]. 



3 



3. Definitions and Notation 

Mathematical definitions, notation and a fundamental supportive theo- 
rem are presented in this section. More specialized definitions will follow in 
appropriate sections. The standard Euclidean norm will be denoted by || ||. 



Definition 3.1. The parameterized Bezier curve, denoted as B(t), of de- 
gree n with control points P m G M 3 is defined by 



Bit) = J2 B m,n(t)Pm,t G [0,1], 



m=0 



n 



where B m ^ n {t) = [J t m {l — t) n 111 and the PL curve given by the points 

{Pq, P\, . . . , P n } is called its control polygon. When Pq = P n , the control 
polygon is closed. Otherwise when P ^ P n , it is open. 

In order to avoid technical considerations and to simplify the exposition, 
the class of Bezier curves considered will be restricted to those where the 
first derivative never vanishes. 

Definition 3.2. A differentiate curve is said to be regular if its first deriva- 
tive never vanishes. 

Definition 3.3. A curve is said to be simple if it is non- self -intersecting. 



The Bezier curve of Definition 3.1 is typically called a single segment 



Bezier curve, while a composite Bezier curve is created by joining two or 
more single segment Bezier curves at their common end points. 

We use B to denote a simple, regular, C 1 , composite Bezier curve 
in IR 3 , throughout the paper. 



Definition 3.4. [21] Let X and Y be two non-empty subsets of a metric 
space (M,d). We define their Hausdorff distance fi(X,Y) by 

u(X,Y) := maxjsup inf d(x, y),sup inf d(x,y)}. 
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Subdivision algorithms are fundamental for Bezier curves [Hj and a brief 
overview is given here. Figure [T] shows the first step of the de Casteljau 
algorithm with an input value of | on a single Bezier curve. For ease of 
exposition, the de Casteljau algorithm with this value of | is assumed, but 
other fractional values can be used with appropriate minor modifications to 
the analyses presented. The initial control polygon P is used as input to 
generate local PL approximations, P 1 and P 2 , as Figure 



la 



shows. Their 

union, P 1 UP 2 , is then a new PL curve whose Hausdorff distance is closer 
to the curve than that of P. 




Subdivision 
> 




P 



ue) and P\ 



reen 



(a) Initial and resultant curves 




(b) Subdivision process 
Figure 1: A subdivision with parameter | 



A summary is that subdivision proceeds by selecting the midpoint of each 
segment of P and these midpoints are connected to create new segments, as 
Figure 1(b) shows. Recursive creation and connection of midpoints continues 
until a final midpoint is selected. The union of the segments from the last 
step then forms a new PL curve. Termination is guaranteed since P has only 
finitely many segments. 

After i iterations, the subdivision process generates 2* PL sub-curves, 
each being a control polygon for part of the original curve j8] , which will be 
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referred to as a sub-control polygon^, denoted by P k for k = 1, 2, 3, . . . , 2\ 
Each P k has n points and their union [J k P k forms a new PL curve that 
converges in Hausdorff distance to approximate the original Bezier curve. 
The Bezier curve defined by [J k P k is exactly the same Bezier curve defined 
by the original control points {P , Pi, . . . , P n } [10J. So [j k P k is a new control 
polygon of the Bezier curve. 

Exterior angles were defined [16] in the context of closed PL curves, but 
are adapted here for both closed and open PL curves. Exterior angles unify 
the concept of total curvature for curves that are PL or differentiable. 
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(a) Orientated lines (b) Vectors 

Figure 2: An exterior angle 



Definition 3.5. The exterior angle between two oriented line segments, 
denoted as P m ~\P m and P m P m+ i, is the angle between the extension oj ' P m -\P m 
and P m P m+ i, as shown in Figure 2(a), Let the measure of the exterior angle 
to be a m satisfying: 

< a m < 7r. 

This definition naturally generalizes to any two vectors, Vi and v 2 , by joining 
these vectors at their initial points, while denoting the measure between them 
as r](vi,V2), as indicated in Figure 2(b). 

Definition 3.6. The curvature of a C 2 curve C(t) parametrized on [a,b] 
is given by 

\C'(t)xC"(t)\\ 



K (t) 



t E [a, b]. 



(1) 



3 Note that by the subdivision process, each sub-control polygon of a simple Bezier 
curve is open. 
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Its total curvature [7] is the integral: dt. 

Total curvature can be defined for both C 2 and PL curves. In both cases, 
the total curvature is denoted by T K (-). The unified terminology is invoked in 
Fenchel's theorem, which follows and is fundamental to the work presented 
here. 

Definition 3.7. [16] The total curvature of a PL curve in IR 3 is the sum 

of the measures of the exterior angles. 

Fenchel's Theorem [7] presented below is applicable both to PL curves 
and to differentiable curves. 

Theorem 3.1. [7J Fenchel's Theorem] The total curvature of any closed 
curve is at least 2tt, with equality holding if and only if the curve is con- 
vex. 

Denote a PL curve with vertices {Pq, Pi, ■ ■ ■ , P n } by P, and the uniform 
parametrization [20] of P over [0, 1] by Z(P)ro ji. That is: 

Z(P) [0 ,i](-) = P, for j = 0,1,.- • ,n 
n 

and Z(P)[o,i] interpolates linearly between vertices. 

Definition 3.8. Discrete derivatives [20] are first defined at the param- 
eters tj = -, where 

Kl'huUj) I) 

for j — 0, 1, ■ ■ • , n — 1. Let 



P> = l'(P) m (t 3 ) 



Pj+i - Pj 



Denote P' = (Pq, P{, ■ • • , P^-x)- Then define the discrete derivative for i(P)[o,i] 



as: 



l'(P) m =l(P') m . 



Intuitively, the first discrete derivatives are similar to the slopes defined 
for univariate real- valued functions within an introductory calculus course. 
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4. Angular Convergence under Subdivision 



Using the notation introduced in Section [3j for a sub-control polygon 
P k , we consider the following analysis, where, for simplicity of notation, we 
repress the superscript and denote this arbitrary sub-control polygon simply 
as P, where P has the corresponding parameters of the indicated control 
points by to,t±, ■ ■ ■ ,t n . And let l(P,i) be the uniform parameterization [20] 
of P on [^i, J] k e {1, 2, 3, . . . , 2'}. That is 

l(P,i) = ^(-P)rfc-i fc.i and l(P,i)(t m ) = P m for m = {0, 1, • • • , n}, 



where t T 



fc-i 

2* 



n2' 



. Note from the domain of l(P,i) that 



t n — tr 



Furthermore, let 



and tm — t 



m—l 



n2 i 



for m — {1, • • • , n}. 



(2) 



ai, a 2 , 



«n-l 



be the corresponding measures of exterior angles of P (Definition 3.5). Treat 
the discrete derivative l'(P,i)(t) at some t G [0,1] as a point in IR 3 , then 
the Euclidian distance between two such discrete derivatives refers to the 
Euclidian distance between two points. Consider the Euclidian distance 
between the discrete derivatives of the two consecutive segments, that is 
\\l'(P,i)(t m ) — l'(P,i)(t m -x)\\. We will show a rate of O(^) for the conver- 
gence 

\\l'{P,i){t m ) - r(P,i)(*m-i)|| -> as i -> oo. 

This will imply that cos(a m ) — > 1 with the same rate and that a m — > at a 
rate of 0( 

Lemma 4.1. For a C 1 , composite Bezier curve B, the value \ \l'(P,i)(t m ) 
l'(P,i)(t m -i)\\ converges uniformly for allm and k to zero at a rate ofO(^ 



Proof: 



\l'(P,i)(t m )-l'(P,i)(t 



m—l) 



< \\l'(P,i)(t m )-B'(t m )\\ + ||£'(i m ) - B'(t 



m—l ) 



+ \\B'(t. 



m—l , 



l'(P,i)(t 



m—l ) 



(3) 
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The first and the third terms converge to at a rate of O(i-), because the 
discrete derivative converges to the derivative of the original curve with this 
rate M- 

For ease of exposition in this proof, we will initially assume that B is 
C 2 , but then complete the proof by extending this special case to the more 
general assumption that B is C 1 . If B is C 2 , then for the second term, the 
Mean Value Theorem on (0, 1) provides that there exists some s G (t m ,t m -i) 
such that 

\\B'{t m )-B'{t m ^)\\ = \\B"(s)\\ ■ \t m -t m ^\. (4) 

Let 7 = supj 0jl ] ||B"(t)||, which exists because of the continuity of B" over 
[0, 1]. It follows that 

\\&{t m )-B , (t m -l)\\<8U V \\& , (t)\\-\t m -t m - 1 \ = r ~r\A ' (5) 



[0,1] 



n2 l 



The second equality holds by Equation (|2]). So | \l'(P, i)(t m ) — l'{P, i)(t m -i)\\ 
converges to zero at a rate of O(^). 

The curve B is only assumed to be C 1 at its junction points, so it is natural 
to question the validity of Equation [4] and Inequality [H] in the C 1 . Note that 
Equation [4] is over the open interval (0,1), so no difficulties occur relative 
to the junction points. However, more subtle consideration is necessary for 
Inequality [5} Each segment of the composite Bezier curve is a polynomial 
defined over [0, 1] , so the definition of 7 as a supremum over [0, 1] still exists. 

Theorem 4.1 (Angular Convergence). ForaC 1 , composite Bezier curve 
B, the exterior angles of the PL curves generated by subdivision converge uni- 
formly to at a rate of 0(< 

Proof: Since B(t) is assumed to be regular and C 1 , the non-zero minimum 
of ||£>'(t)|| over the compact set [0,1] is obtained. For brevity, the notations 
of Ui = l'(P,i)(t m ), Vi = l'(P,i)(t m -i) and a = a m (measure of an arbitrary 
exterior angle) are introduced. The convergence of Ui to B'(t m ) [20J implies 
that \\ui\\ has a positive lower bound for i sufficiently large, denoted by A. 

gives that ||itj — Vi\\ — > as i — > 00 at a rate of 0(4)- This 



Lemma 



4.1 



implies: ||itj|| — \ \vi\\ — > as i — > 00 at a rate of 0(i 

Consider the following where each multiplication between vectors is the 
dot product: 



1 — cos(a) = 1 — r; — -j — j-j 
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- ~\ + II II ^ x + \ ^ \ 6 
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kll 




kl 





It follows from Lemma 4.1 that the right hand side converges to at a rate 
of O(i). Consequently by the above inequality 1 — cos(a) — > with the 
same rate. It follows from the continuity of arc cos that a converges to as 
i — > oo. 

To obtain the convergence rate, taking the power series expansion of cos 
we get 

l-co 8 (a)>(a)Ml-| 1 f- 1 f + -|) 



Note that for 1 > a 



2 



111 1 a 

e=1 + 1+ 2! + 3! + 4! + '" >l iri3r + '"l' (8) 
Combining Inequality [7] and |8] we have, 

1 — cos(a) > (a) 2 ■ (- — a 2 ■ e). 

For any < r < |, sufficiently many subdivisions will guarantee that a is 
small enough such that 1 > a and r > (a) 2 • e. Thus 

1 - cos(a) > (a) 2 • (- - a 2 ■ e) > (a) 2 • - r) > 0. 

Then convergence of the left hand side implies that a converges to at a 
rate of 

5. Topological Equivalence (Homeomorphism) 

The terminology of topological equivalence is classical, referring to two 
sets being homeomorphic. We present sufficient conditions for a homeomor- 
phism between a subdivided control polygon and its associated Bezier curve. 
To establish a homeomorphism under subdivision, we first establish a local 
homeomorphism between a sub-control polygon and the corresponding sub- 
curve of £>, and then establish a global homeomorphism between the control 
polygon and B. 
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5.1. Local arguments for topological equivalence 

A similar lemma has appeared [23] without proof. It is likely that the 
proof is well known, possibly as a 'folk theorem', but we are unable to find 
a reference in the literature, so its proof is provided here, for the sake of 
completeness. 




Pj(Qs) 



>P4 (Qs) 

P~5 (Q«) 

(a) A single point intersection 




P2(Q0) 

(b) Coincident line segments 
Figure 3: Self-intersecting PL curves in 3D 



Lemma 5.1 (Non-self-intersection criteria^]). Let P = (P ,Px, ■ ■ ■ ,P n ) 
be an open PL curve in IR 3 . If T K (P) = Y^jZi a j < 7r ; then P is simple. 



Proof: Assume to the contrary that P is self-intersecting. Then there 
must exist at least one closed loop. (The assumption Y^j=i a j < 71 precludes 
the case of two consecutive edges being coincident.) Arbitrarily choose one 
such loop. There are two cases to consider: a single point intersection as in 
Figure |3(a)| and the coincident case as in Figure 3(b)[ 

Case 1 (A single point intersection): Label the single intersection as Qq 



to be the first vertex of the loop (Figure 3(a)). Label the other vertices 



following the orientation of P and denote this loop by 



Q = (Q ,Q 



i) ■ • • j 



Qn' ) Q 



0), 



for an appropriately chosen value of n'. Denote the measure of the cor- 
responding exterior angle of Q at Qo, Qi, • • • , Q n i by /3q, 0i, ■ ■ ■ , n '- Note 



that [3q < it (Definition 3.5). Note also that each exterior angle (3i (for 
I = 1, 2, • • • , n') of Q (except the initial angle /?o) is an exterior angle of P. 
So 

n' n' n—1 

1=0 1=1 j=l 
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Since YTj=\ a j < 17 we § e ^ 



n-1 



3=1 



1=0 



But for the closed loop Q, XwLoA — 2,71 (Theorem 3.1), which is a contra- 
diction. 



Case 2 (Coincident line segments): The loop is formed by a subset of 
the vertices {P ,Pi,--- ,P n } in this case (Figure 3(b)). Pick an arbitrary- 
point from the subset and relabel it as Qq. Label the other vertices following 
the orientation of P and use the notation Q = (Qa,Qi, ■ ■ ■ ,Q n ',Qo) and 
0o, flx, • • • , j5 n i as the above. Note for this case that, each exterior angle /3/ of 
Q is an exterior angle of P. So 



n-1 



1=0 



3=1 



And YTj=i a j < so S/Lo A < 71 ■ But for the closed loop Q, Y^i=ofii — ^ n 



(Theorem 3.1), which is a contradiction. 



□ 



Theorem 5.1 (Simple sub-control polygons). ForaC 1 , composite Bezier 
curve B, there exists a sufficiently large value ofi, such that after i -many sub- 
divisions, each of the sub-control polygons generated as P k for k = 1, 2, 3, . . . , T 
will be simple. 

Proof: For each P h , the measures of the exterior angles of P k converge 
uniformly to zero as i increases (Theorem 4.1). Each open P k has n edges. 
Denote the n — 1 exterior angles of each P k by a k , for j = 1, . . . , n — 1 and 
for k = 1, 2, 3, . . . , 2\ Then there exists i sufficiently large such that 



n-1 

E 

3=1 



OLj < 7T, 



for each k — 1, 2, 3, . . . , 2\ Use of Lemma 5.1 completes the proof. 



□ 
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5.2. Global arguments for topological equivalence 

The proof techniques for homeomorphism rely upon the sub-control poly- 
gons to be pairwise disjoint, except at their common end points. Denote two 
generated sub-control polygons of B as 

P=(P ,P U --- ,P n ) aadQ = (Q 0> Q 1> ... ,Q n ). 

Definition 5.1. The sub-control polygons P and Q are said to be consec- 
utive if the last vertex P n of P is the first vertex Qo of Q, that is, P n = Qq. 



Remark 5.1. For B, the C 1 assumption ensures that the segments P n ~\P n 
and QoQi are collinear, so the exterior angle at the common point P n = Qq 
is either or tt. But the regularity assumption ensures that the exterior angle 
can not be n, so that the exterior angle at the common point is 0. 



Lemma 5^ extends to arbitrary degree Bezier curves from a previously 
established result that was restricted to cubic Bezier curves |25j, as used in 
the proof of isotopy under subdivision for low-degree Bezier curves 




Figure 4: Intersecting consecutive sub-control polygons 



Lemma 5.2. Let IT be the plane normal to a sub-control polygon at its initial 
vertex. If the total curvature of the sub-control polygon is less than |, then 
the initial vertex is the only single point where the plane intersects the sub- 
control polygon. 
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Proof: Denote the sub-control polygon as Q = (Q , Qi, . . . , Q n ), where 
Figure [4] shows an orthogonal projection of this 3D geometry. Assume to the 
contrary that II fl Q contains a point u where u ^ Qo- Assume without loss 
of generality that a segment Q m -\Qm intersects II at the point u for some 
m G {2, 3, ■ ■ ■ , n}. This would produce a closed PL curve denoted as 

J = (Qo, Ql, • • • , Qm-l,U, Qo)- 

Denote the exterior angle of J at Qq as 9 and note that 9 = |. Denote the 



3.1 we have 



exterior angle of J at u as 7. Then 7 < tx (Definition 3.5). Denote the other 
exterior angles of J as (3i,i = 1, . . . , m — 1. By Theorem 

m—1 

9 + J2 A + 7 > 2tt. 



i=i 

Thus 

m—1 _ 

^/3 i >27r-0-7> \ 
i=i 

These fa are also exterior angles of Q resulting in lower bound for the total 
curvature of Q as given by T K (Q) > Pi > vr/2, which is a contradiction 

to T K {Q) < f . ^ □ 

Lemma 5.3. Recall that B denotes a simple, regular, C 1 , composite Bezier 
curve in M. 3 . Let w be a point of B where B is subdivided and let II be the 
plane normal to B at w. Then there exists a subdivision of B such that the 
sub-control polygon ending at w and the sub-control polygon beginning at w 
intersect IT only at the single point w. 

Proof: The plane II separates M. 3 into two disjoint open half-spaces, 
denoted as Hi and H 2 , such that E 3 = Hi U II U H 2 and Hi n H 2 = 0. 



By Remark 5.1 the exterior angle at {w} is 0. 

Perform sufficient many subdivisions so that the control polygon ending 
at w, denoted by P , and the control polygon beginning at w, denoted by 



Q, each have total curvature less than | by Theorem 4.1 Therefore, by 



Lemma 5.2 the only point where P or Q intersect II is at w. □ 



This global homeomorphism will be proven by reliance upon pipe surfaces, 
which are defined below. 
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Definition 5.2. The pipe surface of radius r of a parameterized curve 
c(t), where t G [0, 1] is given by 



p(t,6) = c(t) + r[cos{9) n(t) + sin(6) b(t)}, 

where 9 e [0, 2tt] and n(t) and b(t) are, respectively, the normal and bi- 
nomial vectors at the point c(t), as given by the Frenet-Serret trihedron. The 
curve c is called a spine curve. 

For B and i subdivisions, with resulting sub-control polygons P k for 
k = 1,. . . ,2\ let Se(r) be a pipe surface of radius r for B so that S r (B) 
is nonsingular. For each k = 1, . . . , 2 l , denote 

• the parameter of the initial point of P k by t^, and that of the terminal 
point by t k n 

• the normal disc of radius r centered at B(t) as D r (t), 

• the union U<e[< fc t k ] D r (t) by T^, and designate it as a pipe section. 

Theorem 5.2 (A simple homeomorphic control polygon). Sufficient 
subdivisions will yield a simple control polygon that is homeomorphic to B. 



Proof: By Theorem 4.1 we can take l% subdivisions so that T K (P k ) < tt/2, 
for each sub-control polygon P k . By Lemma 5.1 this choice of i\ guarantees 



that each P k is simple. By the convergence in Hausdorff distance under 
subdivision |22j, we can take subdivisions such that the control polygon 
generated by i% subdivision fits inside the pipe surface S r (B). Choose i = 



max{ti, 12}. By Lemma 5.3, this choice of 1 ensures that each P k fits inside 
the corresponding Y\.. This plus the fact that P k is simple shows that the 
control polygon, U fc=1 P , is simple, which implies the homeomorphism. □ 



6. Sufficient Subdivision Iterations 

In this section, we shall establish: (1) numbers of subdivisions for small 
exterior angles based on Angular Convergence; and (2) numbers of subdivi- 
sions for homeomorphism. 

From the previous sections we know that the homeomorphism is obtained 
by subdivision based on two criteria: (1) Angular Convergence; and (2) con- 
vergence in distance. So the speed of achieving this topological characteristics 
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is determined by the Angular Convergence rate and the convergence rate in 
distance which are both exponential. Here, we further find closed-form for- 
mulas to compute sufficient numbers of subdivision iterations to achieve these 
properties. 

Definition 6.1. Let P denote a control polygon of a Bezier curve, and let 
P x denote an ordered list of all of x- coordinates of P (with similar meaning 
given to P y for the y- coordinates and to P z for the z- coordinates). Let 

|| ^%Px || oo max \Pm—l,x 2P m:r -\- Pjn-\-\ x\ 
0<m<n 

be the maximum absolute second difference of the x-coordinates of control 
points, (with similar meanings for the y and z coordinates) . Let 

= (|| A 2 -Pr Hoc, || A 2 P y || || A 2 P Z \\oo), 

(i.e.) a vector with 3 values. 

Definition 6.2. The parameter measure distance^ [22] between a Bezier 
curve B and the control polygon l(P,i) generated by i subdivisions is given by 

max \\l(P,i)(t) -B(t)\\. 
te[o,i] 

Lemma 6.1. The parameter measure distance between the Bezier curve and 
its control polygon after ith-round subdivision is bounded by 

ijVooCrOIIAaPH, (9) 

where 

= WW], 

V ; 2n 

Proof: A published lemma [221 Lemma 6.2] proves a similar result re- 
stricted to scalar valued polynomials. We consider coordinate-wise and apply 
this result to the x, y, and z coordinates respectively, so that the distance of 



5 The paper |22j uses 'distance', but we adopt this terminology 'parameter measure 
distance' instead. 



16 



the x-coordinates of the Bezier curve and its control polygon after ith-round 
subdivision is bounded by 

^ri\Un) || A 2 P X \U 

with similar expressions for the y and z coordinates. Taking the Euclidean 
norm of the indicated three x, y and z bounds yields the upper bound given 
by (J9|, an upper bound of the parameter measure distance between the Bezier 
curve and its control polygon after the ith subdivision. □ 



For convenience, denote the above bound in distance as: 



B**t(i) ■= ^oc(n)||A 2 P||. (10) 



2 

Lemma 6.2. After i subdivision iterations, the distance between l'(P,i) and 
B' is bounded by B' dist {i), where 

B' dist (i) := ^rN^n - 1)\\A 2 P% (11) 
and P' that consists of n — 1 control points is the control polygon of B' . 

Proof: The derivative l'(P,i) of the control polygon for a Bezier curve is 
identical to the control polygon l(P', i) for the Bezier curve's first derivative 
by [20~t Lemma 6]. That is, 

l'(P,i) = l(P',i). 



Since (Lemma 6.1) 



we have 



^\\l(P',i)(t)-B'(t)\\<B' dist (i), 
te[o,i] 



m^\\l'(P,i)(t)-B'(t)\\<B' dist (i). (12) 
te[o,i] 

□ 
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6.1. Sufficient subdivision iterations for small exterior angles 

Assume v is a small measure of angle between and 7r. We shall find how 
many subdivisions will generate a control polygon such that the measure a 
of each exterior angle satisfies 



a < v. 



(13) 



Recall the proof of Angular Convergence (Theorem 4.1). Consider two 
arbitrary consecutive derivatives Ui = l'(P, i)(t m ) and Vi = l'(P,i)(t m -i) 
and the corresponding exterior angle a. Recall that in Section [4] we had 
Inequality [6] 

2|k - Ui\\ 



1 — cos(a) < 



U; 



Recall the proof of Lemma 4^ where Inequality [3] is 



Ui - Vi 



< ||r(p, i )(U-sxuil + ll^(U-s / (Wi)ll + l|i3 / (t m -i)-/ / (^0(^-0 



and Inequality [5] is 

| \B\t m ) - B'{t m -i) 1 1 < sup \\B"{t) || • \t m - * m _i| 

[o,i] 



7 

n2 i 



where 7 = supr 01 i Combining the above inequalities yields 

. . 2 
1— cos(a) < 



2(\\l'(P,i)(t m )-B'(t m )\\ + \ 


\B (t m 


_ 1 )-l'{P,i){t m _ l )\\+ 1 /{n2 i )) 




\Ui\ 





Using Lemma 6.2 we get 



1 — cos(a) < 



2{2B' dist {i) + l/{n2 i )) 



Ui 



(14) 



Let a = min{||i3'(t)|| : t G [0,1]}. The regularity of B ensures that 
a > and the continuity of £>' on the compact interval [0, 1] ensures that the 
minimum exists. Recall Ui = l'(P,i)(t m ) for some t m 6 [0,1]. So it follows 



from Inequality 12 that 



\B'(t r 



Ml < B'distd)- 
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Solving the inequality we get 

|kli>iiB'(Uii-^k t «>^-^k^)- 

In order to have ^ 0, it is sufficient to perform enough subdivisions such 
that 

IN I > cr - B' dist (i) > 0, 



that is B' dist (i) < a. By the definition (Equation 11) of B' dist (i) we set, 

1 



2 2 



A 2 P' ||< o. 



Therefore for B' dist {i) < a, it suffices to hav^] 



. 1 Mn-1) || A 2 P' 

i > - log( 

Z a 



(15) 



It is worth noting that Ni is a subdivision number depending on variables 
P', n and a. 

After the i subdivision iterations, whenever i > Aq, then B' dist {i) < 
B dist (Nx), because B' dist (i) is a strictly decreasing function (Equation 11). 
So it follows from Inequality 14 that whenever i > Aq, 

! _ cos(a) < 2(2B^(i)_+ 7 /(n2')) < 2(23^(0 + 7 /(n2<)) 



To obtain a < v (Inequality 13), it suffices to have that 1 — cos(a) < 
1 — cos{y). Now choose i large enough so that 

° - B dist\ N V 

The second inequality of Inequality [16] implies that 



2B' dist (i) + -L < 1(1 - cos(i/))(a - B'^N,)). 

We could solve the above inequality for i, but we can avoid this compli- 
cated computation by noting that 2B' dist {{) is much smaller than when 



Throughout this paper, we use log for log. 
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i is large. (Note that 7 = supr 01 i ||jB"(i)||. If 7 = 0, then the Bezier curve 
would be a straight line segment. We exclude this trivial case and assume 
that 7 > 0.) 

So let 2B' dist (i) < ^, that is 

2nN^n-l) || A 2 P' || = 

7 

So for i > N2 we get 

diStV ; n2* n2 J 



Then we consider 



2^- < ^(1 - cos{v)){a - B^W). (18) 



Solve Inequality 18 and get 



1 Tl 

-<-(l-cos(z,))(a- J B^(iV 1 )). 



So 

2 > lo. 



na-cos^))^-^^))^ 
To simplify this expression, let 

47 

CB= n{a-B> dM {N 1 )y 
which is a constant for a given £>, and let 

fiy) = r^^- (is) 

1 — COS V 

Then, we have 

% > log(/(i/)). 

Theorem 6.1. Given any v > 0, there exists an integer N{u) defined by 

N{u) = maxima, N 2 , Iog(/(z/))} (20) 



where N±, N 2 and f{v) are given by Equations \15 11 and 19 respectively, 
such that each exterior angle is less than v, whenever i > N(u). 

Proof: It follows from the definitions of N\, N 2 and f{y) and the analysis 
in this section. □ 
It is worth to note that N is a logarithm depending on several parameters 
such as 7, a, iVoo(n) and A 2 P' as well as an upper bound variable v. 
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6.2. Sufficient subdivision iterations for homeomorphisms 

For a Bezier curve B of degree 1 or 2, the control polygon is trivially simple 
and ambient isotopic to B, provided the regularity of B. So we consider n > 3. 



Given any v > 0, Theorem 6.1 shows that there exists an integer N(u), 
such that each exterior angle is less than v after N(u) subdivisions. Further- 
more, there is an explicit closed formula to compute N(u). 



There exists a positive integer, iV(^-j- 
such that after N(^- 



Theorem 6.2 

iV(^-j-) is defined by Equation 
sub-control polygon will be simple. 
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for n > 2, where 
) subdivisions, each 



Proof: By Theorem 
less than 



6.1 



after N(-^j) subdivisions, each exterior angle is 
. Since each sub-control polygon has a n — 1 exterior angles, the 
total curvature of each sub-control polygon is less than tt. Lemma [oTT] implies 
that this is a sufficient condition for each sub-control polygon being simple. □ 



While existence of sufficiently many iterations for the control polygon to 
fit inside the pipe S r (B) has been established, it remains of interest to bound 
the number of subdivisions that are sufficient for this containment. Define 
JV'(r) by 



1 



N'(r) = - log 



.AU™)I|A 2 P| 



(21) 



2 r 

where r is the radius of a nonsingular pipe surface for B. By the definition 
of B<nst{i) (Equation 10) and Equation 21, we have B^^i) < r whenever 
z > N'(r). 

Lemma 6.3. The control polygon generated by i subdivisions, where i > 



N'(r) and N'(r) is given by Equation 21, satisfies 



max \\B(t) -l(P,i)(t)\\ < r, 

te[o,i] 

and hence fits inside the pipe surface of radius r for B. 



Proof: By Lemma 6.1, maxt e [o,i] \\B(t) — l(P,i)(t)\\ < Bdi S t(i)- Then this 



lemma follows from the definition of N'(r) given by Equation 21 □ 
While Theorem £^2 addresses each sub-control polygon, it is of interest 
to ensure that the union of all these sub-control polygons is also simple. In 
Theorem |6.3[ that union is the 'control polygon', as the result of multiple 



subdivisions. 



21 



Theorem 6.3. Take any N > max{iV( 



2(n-l). 



fined by Equations 20 and N'{r) is given by Equation 21 
subdivision, the control polygon will be simple. 



,N'(r)}, where N(u) is de- 
Then after the Nth 



Proof: The inequality iV > N'(r) implies that the control polygon 
generated after the Nth subdivision lies inside the pipe. The inequality 
N > -M"( 2 ( n -i) ) ensures that the total curvature of its each sub-control poly- 
gon is less than |. These two conditions are sufficient conditions for the 
control polygon being simple (The proof of Theorem 5.2). □ 



7. Conclusion 

We proved the Angular Convergence and used it to show homeomorphism 
between a Bezier curve and its control polygon under subdivision. We showed 
that the rate of Angular Convergence is exponential. Plus the exponential 
rate of convergence in terms of Hausdorff distance, we established closed-form 
formulas to compute sufficient numbers of subdivision iterations to achieve 
topological properties. 

The Angular Convergence is shown to be useful for determining homeo- 
morphism. It is natural to consider ambient isotopy next [T2] , where total 
curvature will be fundamental, according to Fary-Milnor Theorem [T6] . 
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